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' Abstract. Meridional flow results from slight deviations from the thermal wind balance. The 
(n: deviations are relatively large in the boundary layers near the top and bottom of the convection 
zone. Accordingly, the meridional flow attains its largest velocities at the boundaries and de- 
creases inside the convection zone. The thickness of the boundary layers, where meridional flow 
is concentrated, decreases with rotation rate, so that an advection-dominated regime of dynamos 
' is not probable in rapidly rotating stars. Angular momentum transport by convection and by the 
^SJ . meridional flow produce differential rotation. The convective fluxes of angular momentum point 
radially inward in the case of slow rotation but change their direction to equatorward and paral- 
' , lei to the rotation axis as the rotation rate increases. The differential rotation of main-sequence 
Pm . dwarfs is predicted to vary mildly with rotation rate but increase strongly with stellar surface 
^0 ■ temperature. The signiflcance of differential rotation for dynamos has the opposite tendency to 
(-H ] increase with spectral type. 
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. 1. Introduction 

Differential rotation and meridional flow are the two main components of global stel- 

1- H ■ lar circulation. Both are important for dynamos. Differential rotation winds toroidal 
K>" I magnetic fields, and meridional flow near the base of the convection zone is probably 

responsible for the observed equatorward migration of sunspot activity. Knowledge of 
differential rotation as a function of stellar parameters is a key for understanding stellar 
■ dynamos. 

, Helioseismology has shown that the regions inside the sun occupied by thermal convec- 

I ■ tion and differential rotation coincide. The same is probably true of the meridional flow. 
' Any theory of differential rotation and meridional flow has to describe the global flows 
T-H ■ against the background of convective turbulence. The large-scale and turbulent flows are 
^ \ intimately linked. The theory, therefore, has to rely on the tools and methods of the 
•j^ • mean-field hydrodynamics of turbulent fluids. 

Kn] ] This paper reviews the mean-field theory of meridional flow and differential rotation. 

The differential rotation can be understood as a result of interaction between convection 
and rotation. The meridional flow is produced by the non-conservative part of centrifu- 
gal force and by the baroclinic torque provided that these two drivers do not balance 
each other. Numerical models based on the theory are also discussed. The models repro- 
duce closely solar rotation and predict the dependence of differential rotation on stellar 
parameters. The results of 3D numerical simulations are discussed only occasionally. Nu- 
merical experiments have been reviewed recently by Miesch & Toomre (2009) and Brun 
& Rempel (2009). 
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2. Meridional circulation 

Meridional circulation is a vortical flow. Accordingly, the flow is convenient to describe 
in terms of the azimuthal vorticity uj = (V x V")^, where V™ is the global (azimuthally- 
averaged) meridional velocity. The mean-fleld vorticity equation, 

'^+rsineV.(v-^)+ViVn- rf! - , (2.1) 

ot \ rsmO/ oz CpV of) 

collects the meridional flow drivers on its right side while the left side describes the 
flow reaction to this driving. In particular, the term V (V") accounts for the viscous 
resistance by the eddy viscosity to the meridional flow driving, 

ViV-') (-_(,M,.^)j, (2^2) 

where Afkim is the eddy viscosity tensor and repetition of subscripts signifies summation. 
In Eq. (|2.ip . the usual spherical coordinates (r, 9, (f)) are used, is the angular velocity, 
S is the specific entropy, g is gravity, Cp is the specific heat at constant pressure and 
d/dz = cos9d/dr — sinOd/dO is the spatial derivative along the rotation axis. 

In hydrodynamics, there are only two main sources of meridional flow illustrated by 
Fig.[TJ The first term on the right side of Eq. (12. ip shows that differential rotation can 
produce the flow (Kippenhahn 1963). If angular velocity varies with the cylindrical co- 
ordinate z along the rotation axis, the centrifugal force is not conservative and generates 
vorticity. This effect was recently named 'gyroscopic pumping' (Garaud & Bodenheimer 
2010). If angular velocity decreases with distance from the equatorial plane, as it does 
in the sun, the centrifugal force produces a torque driving anti-clockwise circulation (in 
the north-west quadrant of the meridional cross-section) . 




Figure 1. Illustration of centrifugal {left) and baroclinic {right) driving of the meridional flow 

(see text). 

The second term on the right side of Eq. (|2.ip allows for the baroclinic driving of the 
meridional flow. If polar regions are warmer than the equator, as they probably are on 
the sun (Rast et al. 2008), then the warm fluid tends to rise up near poles and spread 
over the surface while the cooler equatorial fluid tends to sink down and spread over near 
the base of the convection zone thus producing clockwise circulation (Fig.[T]). 

The two sources of the meridional flow are competing on the sun. This is probably 
not by chance. It seems to be a general rule that these two sources are nearly bal- 
ancing each other in convective stars. If Eq. (|2.ip is normalized to dimensionless units 
by multiplying this equation by the square of the viscous diffusion time, R'^/iy'^, the 
first term on the right side gets a coefficient of the Taylor number, Ta = Ail^W^ /v'^; 
R is the stellar radius and is the eddy viscosity. This number is large on the sun. 
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Ta ^ 10^. The second (baroclinic) term also gets a big coefficient of the Grashof number, 
Gr = {gR^/i/^){ST/T) ^ lO'^STsmf, ST is the pole-equator temperature difference and 
STsuri is the surface value of this differential temperature. The left side of the equation, 
on the contrary, scales with the much smaller Reynolds number Re = VqR/v^ ~ 10, 
where Vb ~ 10 m/s is the characteristic amplitude of the meridional flow (the second 
term on the left side of (I2.ip scales with Re^). Therefore, the Eq. (|2.ip has to be satisfied 
by the two terms on the right side nearly balancing each other. 

Neglecting the left side in Eq. (12.11) gives the famous equation of thermal wind balance 

' = ""'^^ - —rW 

The equation shows in particular that a positive differential temperature of the order of 
1 K is necessary for isorotational surfaces to deviate considerably from cylinders in solar 
rotation models. 

Equation (j2.3l) shows what the thermal wind balance is but it does not explain how the 
balance is maintained. In order to understand this, we have to return to the meridional 
flow equation p.l|) and imagine that the balance is somehow violated. This gives a strong 
source of meridional flow. The excited flow reacts back on the distributions of angular 
velocity and entropy by transporting angular momentum and heat to reestablish the 
balance. Meridional flow maintains the thermal wind balance and deviations from the 
balance drive the flow. 
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Figure 2. Lett: depth profile of the meridional velocity for 45° latitude after the mean- field 
model of Kitchatinov & Olemskoy (2011). Negative velocity means poleward flow. Right: depen- 
dencies of the baroclinic (dashed) and centrifugal (dashed-dotted) terms of Eq. (|2.ip and their 
sum (full line) on depth for the same latitude. The meridional flow attains its largest velocities 
in the boundary layers where the thermal wind balance is violated. 

This qualitative picture can explain why the meridional flow decreases with depth 
beneath the solar surface (Zhao & Kosovichev 2004). It does so because the stress- free 
conditions on the boundaries of the convection zone are not compatible with thermal 
wind balance. Accordingly, the boundary layers form where the balance is violated (Dur- 
ney 1989). As deviations from the balance produce the meridional flow, the meridional 
velocity attains its maximum values at the top and bottom boundaries and decreases in- 
side the convection zone. Fig. [5] illustrates the thermal wind balance and the meridional 
flow structure typical of mean-field models of global solar circulation. Boundary layers 
were found also in the 3D simulations of Brown et al. (2008). The bottom flow of Fig. [5] is 
not small compared to the surface flow. The flow, however, decreases rapidly with depth 
beneath the base of the convection zone (Oilman & Miesch 2004). 
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Fig.[2]shows the boundary layers where the thermal wind balance is violated. The bulk 
of the convection zone is, nevertheless, very close to the balance. Even there, however, 
the balance condition is satisfied only on average (Fig. [5] results from the mean-field 
model). Differential rotation is produced by turbulent convection, so that fluctuating 
deviations from the thermal wind balance are unavoidable (Brun et al. 2010). As the 
deviations produce the meridional flow, the flow is expected to fiuctuate considerably with 
time. This qualitative picture was reproduced quantitatively by Rempel (2005) whose 
global circulation model allowed for random fluctuations in the A-effect (cf. the next 
Section for the A-effect definition). The fluctuating A-effect produces small fluctuations 
in differential rotation, which in turn produce much larger irregular variations in the 
meridional flow. The relative amplitude of fluctuations of the meridional flow in Rempel's 
model was about two orders of magnitude larger compared to fluctuations in angular 
velocity. Helioseismologically detected meridional flow indeed shows considerable changes 
from year to year (Zhao & Kosovichev 2004; Gonzales Hernandez et al. 2006) contrasting 
the slight variations in the rotation rate. 

3. Differential rotation 

Helioseismology shows that differential rotation and thermal convection occupy the 
same region inside the sun. The entire convection zone is rotating differentially while 
the rotation inhomogeneity decreases rapidly with depth beneath the convection zone 
(Wilson et al. 1997; Schou et al. 1998). This supports the theoretical concept pioneered 
by Lebedinskii (1941) that the differential rotation results from interaction between con- 
vection and rotation. Convection in a rotating fluid is disturbed by the Coriolis force. 
The back reaction disturbs rotation to make it non-uniform. 

The mean-fleld equation for angular velocity. 



shows that the angular velocity distribution can be modified by angular momentum 
transport, by convection (u) and by meridional flow (V°). 



The ability of convective motions to transport angular momentum even in the case of 
uniform rotation is called the 'A-effect' (Riidiger 1989). The uniformity of rotation is 
mentioned because the eddy viscosity can transport angular momentum also if rotation 
is not homogeneous. The viscosity, however, tends to diminish the differential rotation, 
not to produce it, in contrast with the non-diffusive A-effect. 

Eq. (|3.1|) shows that the azimuthal and meridional convective velocities have to be 
correlated for the A-effect to emerge. The finite cross-correlation requires the convective 
turbulence to possess a preferred direction of anisotropy or inhomogeneity (cf. Kitchati- 
nov (2011) for pictorial discussion of the origin of the A-effect). The radial preferred 
direction is imposed by gravity. The A-effect is allowed for by the non-diffusive part Qfj 
of the velocity correlation tensor Qij = (uiUj). The radial (Q^^) and meridional (Qg^) 
non-diffusive fiuxes of angular momentum for radially-stratified fluid read 




(3.1) 



3.1. The A-effect 




(3.2) 
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where £ is the correlation length, 



Tt^g dS 
15cp dr 



(3.3) 



is the background eddy viscosity, and r is the convective turnover time. The dimensionless 
functions V{il*) and in Eq. (13.2^ depend on the Coriolis number 

(3.4) 



which is the key parameter of the differential rotation theory. The parameter measures 
intensity of interaction between convection and rotation. Its value defines whether the 
convective eddies are living long enough for rotation to influence them considerably. The 
Coriolis number (|3.4p depends on depth in the convection zone. It is smaller than one near 
the solar surface but increases with depth to exceed ten near the base of the convection 
zone. Therefore, the A-efTect theory should be nonlinear in f2*. 

SLOW ROTATION FAST ROTATION 
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Figure 3. The functions V{il*) and H{Q,*) of the A-effect of Eq. p.2p derived by Kitchatinov 
& Riidiger (1993, 2005). The right panel shows the directions of angular momentum transport 
for the cases of slow {fl* <C 1) and fast (f2* ^ 1) rotation. 



The function V(fl*) in Eq. p.2[) is the normalized vertical flux of angular momentum. 
Positive V means radial inward transport. The function H{Q*) is the normalized flux 
parallel to the rotation axis. Its positive value means that the angular momentum is 
transported towards the equatorial plane. The two (not mutually orthogonal) components 
of the angular momentum fluxes are shown in Fig. |3]as functions of the Coriolis number. 
In the case of slow rotation, il* ^ 1, the i7- function is relatively small and angular 
momentum is transported downward in radius. In the opposite case of rapid rotation, 
Q* ^ 1, the y-function is small, and the angular momentum flux is parallel to the 
rotation axis pointing to the equatorial plane. As the Coriolis number increases from 
small to large values, the direction of the non-diffusive flux of angular momentum varies 
smoothly from radial inward to equatorward and parallel to the rotation axis. 

3.2. Eddy viscosities 

After the signiflcance of the correlation of azimuthal and meridional convective velocities 
for differential rotation was recognized, there were attempts to observe this correlation 
using sunspots as tracers (Ward 1965). The flrst measurements gave finite correlation, 
Qg^ cos ^ lO'^m^/s^, with positive Qg^ in the northern and negative in the southern 
hemisphere, indicating the equatorward transport of angular momentum. Subsequent 
measurements of Ribes (1986) and Nesme-Ribes et al. (1993), however, did not confirm 
this result. 
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The equatorial acceleration does not demand a certain sign of Qe^ (a positive value 
of Qe^ cos 9 would mean that angular momentum is permanently pumped towards the 
equator implying an increase of differential rotation with time). The total correlation, 

Q., = 0^+Qr„ (3.5) 

includes the contribution Q'^ of the eddy viscosities in line with the A-effect. The viscous 
fluxes of angular momentum, 

do. ( dn dn\ 

Qrd, = -virsin9— iy2Sm9cos9 rcos9— sm9— , 

^ or \ or o9 J 

= -1^1 sin 61— + sin^ 6* f r cos 6*— - sin 61— j (3.6) 

(Kitchatinov et al. 1994), include two viscosity coefficients, v\ = i/^/i(ri*) and V2 = 
^^t/2(^*)- Rotating convection is anisotropic and the eddy viscosity i>i acting in the 
direction normal to the rotation axis differs from the viscosity vi + i>2 for the direction 
parallel to the axis. 

A positive non-diffusive part, Q^gC0s9 > 0, is required for maintaining equatorial 
acceleration. It is hardly possible, however, to separate the contributions of the A-effect 
from the contribution of the eddy viscosities in the total correlation (13.51) supplied by 
observations. The same separation problem is met by 3D numerical simulations of the 
A-effect. It was, nevertheless, possible to extract the from the numerical simulations 
by Kapyla & Brandenburg (2008) and Kapyla et al. (2011). The analytical results in 
Fig. [3] generally agree with the simulations of the A-effect. 

3.3. Surface shear layer 

The natural boundary condition for the differential rotation problem is the requirement 
that the surface density of external forces be zero, = Qre — 0. The stress-free 
condition means that the global flows are produced by internal processes in the convection 
zone, not imposed externally. The condition on the upper boundary is simplifled by the 
small value of the Coriolis number so that the angular momentum fluxes (13.21) and p.6p 
can be taken in the limit of f2* — > 0. In this limit, the viscosity 1/2 in Eq. p.6|) vanishes 
and vi = v^. The top boundary condition, = 0, then reads 

= - [j^J Vm at r^R. (3.7) 

With positive V^(0) of Fig. [31 Eq. p. 71) requires the angular velocity to increase with 
depth near the surface in at least qualitative agreement with the surface shear layer 
detected by helioseismology (Schou et al. 1998). 

We have seen in Sect. [2] that the meridional flow also has a near-surface boundary 
layer. The near-surface rotational shear and meridional flow are mutually related. This 
relation was recently analyzed by Miesch & Hindman (2011). 



4. Differential temperature 

It is not possible to reproduce the helioseismological rotation law by differential rota- 
tion models if the differential temperature is not allowed for (Brandenburg et al. 1990; 
Riidiger et al. 2005; Miesch et al. 2006). 

Similar to the differential rotation, the differential temperature is currently understood 
as a result of rotational influence on convection. The tensor Xij of the eddy thermal 
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diffusion that controls the convective heat flux, 

^-pTx.o-Q^^. (4.1) 
includes the rotationally-induced anisotropy and quenching 

x„^x,(*(o-w«+0,(o-)eia), (4.2) 

The diffusivity quenching functions of the Coriolis number (13.41) are shown in Fig. 2] 




0.0 h 
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Figure 4. Quenching functions of the eddy thermal diffusivity of Eq. (|4.2|l after the 
quasilinear theory of turbulent transport in rotating fluids (Kitchatinov et al. 1994). 

Directions of the entropy gradient and heat flux (|4.ip do not coincide due to the 
anisotropy of heat transport. Even if entropy varies mainly with radius, the meridional 
flux, 

= pTx,0|| (n*) cos sine—, (4.3) 

does not vanish. For the negative radial gradient of entropy in convection zones, the 
meridional flux ()4.3|) transports heat to the poles. There have been multiple attempts 
to observe differential temperature on the sun. Recent observations of Rust et al. (2008) 
suggest that polar regions are warmer than the equator by about 2.5 K. After the Eq. (I2.3P 
of thermal wind balance, this small differential temperature suffices for the isorotational 
surfaces to deviate considerably from cylinders. 



5. Differential rotation models 

What is usually meant by the 'differential rotation models' are numerical solvers which 
define not only the angular velocity but also meridional flow and entropy distributions 
in stellar convection zones. 

5.1. The Sun 

The solar rotation model in Fig. [5] is close to the helioseismological rotation law. This 
Figure was computed by using two different models - one for the convection zone, and 
another model for the tachocline and radiation zone rotation. Physical conditions in 
convection and radiation zones differ so much that it is not possible to cover both in one 
model. The tachocline modeling did not infiuence the computation of differential rotation 
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in the convection zone in any way but just used the resuhs of this computation as the 
top boundary condition. 




Figure 5. Angular velocity isolines (left) and depth-dependencies of rotation frequencies for 
several latitudes [tight) after the mean-filed model of Kitchatinov & Olemskoy (2011). 

The model of the convection zone rotation in Fig. [S] was tuned with one free pa- 
rameter. Mean-field models can avoid arbitrary prescriptions of any parameters. All the 
modeling needs, including the A-effect, the eddy viscosities and thermal conductivities, 
has been derived using the same quasi-linear approximation of the mean-field theory. 
These derivations, formally, leave no freedom for model design. However, the quasi-linear 
approximation is not expected to be very precise. This raises the question of how sen- 
sitive the model is to variations of basic parameters. The sensitivity is normally quite 
low, except for one parameter. E.g., if the A-effect is reduced by 50% by multiplying 
the non-diffusive fluxes of Eq. p.2[) by a factor of 0.5, the surface differential rotation 
in the solar model decreases by only 10%: from 30% to 27%. It is remarkable that in- 
creasing the A-effect by 50% (factor 1.5) also reduces the differential rotation by same 
10%. Sensitivity to an increase of eddy viscosities or heat conductivities is also rather 
low (decreasing the eddy diffusion can make the model unstable to thermal convection 
and, therefore, inconsistent). The model, however, is rather sensitive to the anisotropy 
of thermal conductivity. If the parameter is introduced in Eq. (|4.2p , 

x^J = X. (^<pms,,+c^^\im^^ , (5.1) 

the resulting differential rotation reacts quite sensitively to variations of this parameter 
and generally increases with C^. Fig. [5] was obtained with = 1.5. This implies that 
the quasi-linear theory underestimates the anisotropy of thermal eddy conductivity. All 
the stellar models discussed below were obtained with the same value of = 1.5. 
Meridional flow for the same model as Fig. [5] is shown in Fig. ^ 

5.2. Solar twins 

Differential rotation has already been measured for many solar-type stars (Korhonen 
2012) motivating stellar applications of differential rotation models. Figures |6] and [7] 
show the differential rotation and meridional flow computed for two stars similar to the 
sun in mass but rotating faster (Kitchatinov & Olemskoy 2011). 

CroU et al. (2006) measured the differential rotation of a solar twin e Eridani using very 
precise photometry of the MOST mission. The star is younger than the sun and rotates 
with a shorter period of Prot — 11 days. The differential rotation and meridional flow 
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Figure 6. Differential rotation and meridional flow computed for solar-type star eEridani. From 
left to right: angular velocity isolines, surface profile of rotation rate, and depth profile of the 
meridional flow at 45°-latitude. 

computed for this star is shown in Fig. [51 The computations give the relative magnitude, 

aDR = l-%^, (5.2) 



of the surface differential rotation of e Eri of odr — 0.13 close to the value of aDR = 0.11 
measured by Croll et al. (2006). The absolute value of the surface differential rotation is 
close to that of the sun. 
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Figure 7. Differential rotation and meridional flow computed for the young rapidly rotating 
dwarf AB Doradus. From left to right: angular velocity isolines, surface profile of rotation rate, 
and depth profile of the meridional flow at 45°-latitude. 

Fig. [7] shows differential rotation and meridional flow computed for the very rapidly 
rotating (Prot = 0.51 day) young dwarf AB Doradus. The computed value of Q!dr = 
4.4 X 10^^ is again close to the differential rotation measured by Donati & Collier Cameron 
(1997) by the method of Doppler imaging. Also in this case the differential rotation in 
dimensional units is close to the solar value. 

The surface profiles of Figs. [5] and [7] are not as smooth as it is on the sun. There 
are peculiarities positioned around the latitude, where the isoline of angular velocity 
tangential to the inner boundary at the equator arrives on the surface. 

The meridional flows of these Figures differ from the flow computed for the sun (Fig. ^ 
by their higher concentration at the boundaries of the convection zone. The boundary 
layers discussed in Sect. 2 become thinner as rotation rate increases. The flow of Fig. [7] 
in the very rapidly rotating star consists of two near-boundary jets linked by very slow 
circulation in the bulk of the convection zone. The meridional flow is believed impor- 
tant for the solar dynamo (Choudhuri 2008). It is, however, hard to imagine that the 
surface-jets flow of Fig. [7] can be significant for dynamos. As a star ages and its rotation 
rate decreases, the boundary layers broaden (Fig. [6]) and the meridional flow eventually 
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transforms to a distributed flow as in Fig. [51 It may be expected that the dynamo of a star 
changes to the advection-dominated type from some other dynamo regime, where merid- 
ional flow is not significant, as rotation rate decreases. This may be the reason for the 
presence of two separate branches for slow and fast rotators in the observed dependence 
of the stellar activity cycle periods on rotation rate (Saar & Brandenburg 1999). 

5.3. Mass and rotation rate dependence 

Applications of the mean-field model to individual stars do not always agree with ob- 
servations so closely as in the two examples given above. However, no one case of clear 
disagreement has been met so far. This motivates using the model to explore the depen- 
dence of differential rotation on stellar parameters. Fig. [8] shows the surface differential 
rotation computed for main-sequence stars of different mass and age. The rotation peri- 
ods of this Figure were specified using Gyrochronology (Barnes 2007, 2009). 




0.00 k , , , , A 0.0 F. , , ^ 
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P„i (days) P^„, (days) 

Figure 8. Surface differential rotation as a function of the rotation period for main-sequence 
dwarfs of subsolar (left) and supersolar (right) masses (from Kitchatinov & Olemskoy 2012). 

Fig.[5]suggests that the differential rotation of a star of given mass varies by about 30% 
only as the star ages. A much stronger tendency is the increase in differential rotation 
with stellar mass. The same tendency for rapidly rotating stars was observed by Barnes 
et al. (2005). 




1 10 

Pro, (days) 

Figure 9. Isolines of the Cn dynamo number (|5.3|l from the same computations as Fig. |8] The 
dynamo number increases with rotation rate but decreases with temperature. 

Differential rotation takes part in stellar dynamos by winding the toroidal field lines. 
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The computations in Fig. |S] predict that the largest differential rotation belongs to the 
hottest F-stars in this Figure. The question arises whether the strong differential rotation 
implies high dynamo activity. The computations suggest a negative answer. The efficiency 
of differential rotation in generating magnetic fields can be estimated by the modified 
magnetic Reynolds number that in dynamo theory is conventionally notated as Cq, 

On = («) 

where Ail, is the angular velocity variation within the convection zone, H is the convection 
zone thickness, and 77^ is the magnetic eddy diffusivity. This number estimates the ratio 
of the rate Ail of magnetic field production by differential rotation to the rate /r}^ of 
diffusive decay of the field. Fig. IHlshows isolines of Cn on the plane of effective temperature 
and rotation period, obtained in the same computations as Fig. |51 Cji increases with 
rotation rate. This is mainly due to the rotational quenching of the eddy diffusivity. The 
dynamo number (j5.3l) increases also with decreasing temperature. This is mainly because 
convection is slower and r/^ is smaller in cooler stars. The increase of 77^ with increasing 
temperature overpowers the increase in differential rotation, so that the dynamo-number 
(|5.3|) decreases. 

The large rotational shear of F-stars is less efficient in winding magnetic fields than 
the small differential rotation of K-dwarfs. The widely spread opinion that low mass 
stars host a^-dynamos producing non-axisymmetric global fields is not supported by 
the estimations of the Co dynamo-number. The small but dynamo-efficient differential 
rotation of these stars favors axisymmetric global fields. This might be the reason why 
the magnetic structure of the M-star observed by Donati et al. (2006) was close to axial 
symmetry. 



6. Summary 

We have seen that meridional flow results from (slight) deviations from the thermal 
wind balance. The balance, in turn, is maintained by the flow. Meridional flow attains 
its largest velocities in the boundary layers near the top and bottom of the convection 
zone where deviations from the thermal wind balance are relatively large. The thickness 
of these layers decreases with increasing rotation rate. The layers in rapidly rotating 
stars are so thin that an advection-dominated regime of the dynamo in these stars is not 
probable. 

Differential rotation is produced mainly by convection and also by the meridional flow 
transporting angular momentum. Differential temperature with poles warmer than the 
equator is very significant for differential rotation formation. Differential temperature 
is currently understood as a result of rotationally induced anisotropy of convective heat 
transport. The anisotropy is the parameter to which differential rotation models are most 
sensitive. 

Differential rotation models reproduce closely the helioseismological rotation law and 
the observed differential rotation of several individual stars. The models predict the 
differential rotation of a star of given mass to vary mildly as the star ages and its rotation 
rate decreases. Differential rotation is, however, predicted to increase strongly with stellar 
mass. The significance of rotational shear for dynamos estimated with the Co dynamo 
number (j5.3p has the opposite tendency. The large differential rotation in F-stars is less 
dynamo-efficient than the tiny rotation inhomogeneity of M-stars. 
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